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Abstract 

We perform ab initio electronic structure calculations of the intermetal- 
lic compound CeFeGe3 by means of the Tight Binding Linear Muffin-Tin 
Orbitals-Atomic Sphere Approximation (TB-LMTO-ASA) within the Local 
Spin Density Approximation containing the so-called Hubbard correction term 
(LSDA+U S/C ), using the Stuttgart's TB (Tight Binding)-LMTO-ASA code in 
the framework of the Density Functional Theory (DFT) . 
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1 Introduction 



The heavy electron compounds refer to those having a specific heat coefficient 
7 of the order of Jmol _1 K -2 , which is much larger than that of simple me- 
tals (being typically in the range of mJmol _1 K -2 ). In addition, this implies 
that such compounds— which have been given a great deal of attention for long 
time— possess a large effective mass m*, outweighing several hundreds of times 
the free electron mass [1]. They can be classified into two groups: the concen- 
trated Kondo compounds CK and the intermediate valence compounds (IV) [2,3] 
depending on the position of the 4/ or 5/ level, relative to the Fermi level. The 
CK compounds have an integer valence at temperatures much higher (T > Tt) 
than the Kondo temperature Tfc, at which, there appears the Kondo effect (an 
effect observed in metals with a magnetic impurity). However, at comparatively 
low temperatures (T <^Tk) they form a Fermi liquid state [4] with a reduction 
of its magnetic moment. On the other hand, the IV compounds do not possess 
an integer valence at room temperature as a result of the strong hybridization 
between the 4/ electrons and the conduction electrons, due to the anomalous 
proximity between the 4/ and the Fermi levels. 

Remarkably, the CeFeGe3 compound here studied [5] apparently gives rise to 
two behaviors: A high T& of the order of 100 K and an integer valence for Ce at 
room temperature. This study is motivated by the fact that such compound is 
known to present a strong electronic correlation character, which occurs when 
the Coulomb repulsion between electrons strongly inhibits their motion, thus be- 
coming highly localized. Because of this, on compounds containing lanthanidcs 
the normally expected metal behavior— cerium ions with 4/ electrons or ura- 
nium and neptunium ions with 5/ electrons— is not observed. Some examples 
are: CeAi3, CeCuSi2, CeCu6, UBei3, UCdn, V>-iLx\\-j and NpBei3, a s well as 
transition metal oxides, organic metals and carbon compounds, i.e., carbon 
nanotubes, etc. 

In the first stage of the present work ah initio calculations were carried out 
to investigate the electronic structure of the intermetallic compound CeFeGc 3 
(tetragonal structure with spatial group 107) by using a DFT method [7] in 
the LMTO-ASA approximation [8], whereas in the second, an LSDA+U S/C ap- 
proximation [9] was used, as implemented in the Stuttgart's TB-LMTO-ASA 
code version 47 [8]. The density of states (DOS), total and partial, for cerium 
and iron, as well as the band structure (BS) are obtained from the compound 
geometrical conformation optimized by the CASTEP program (Cambridge Se- 
rial Total Energy Package) which makes use of the ultrasoft pseudopotentials 
theory [10]. The Coulomb and exchange parameters used in the calculations 
were U = 5.4 cV [11] for cerium and U = 2.3 cV [12] (J = 0.9 eV) for iron, just 
obtained in the literature. 
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2 Computational approach. 



2.1 LMTO-ASA and CASTEP Calculations. 



The LMTO-ASA approximation employs the unitary cell splitting into over- 
laping Wigner-Seitz (WS) spheres with a maximum overlap of 15%, which is 
generally considered a reasonable approximation when having a spherically sym- 
metric potential within the spheres. In addition, use is made of the ASA (Atomic 
Sphere Approximation) condition, i.e., a spherical approximation containing no 
zone of free electrons on the Muffin-Tin structure. For open structures, as the 
ones here considered, a set of empty spheres is introduced as a device for des- 
cribing the repulsion potentials in the interstices (between an atomic and an 
interstitial sphere an overlap of 20% is allowed). The total volume of the WS 
spheres is equal to the unit cell volume, thus eliminating the interstitial region. 

On the other hand, the Muffin-Tin (MT) orbital is energy dependent with the 
following linearized form: 

6 (r)-i l Y, (r)f M E ,r)+PA(r/S R ) 1 ; r < S R \ ( 
0aW - i Y lm (r) | {s/r)l+1 ; r>S R j W 

where </> A (r) is found by numerical solution of the radial Schrodinger equation, 
A = RL. R is the site index, whereas I and m are the orbital and magnetic 
quantum numbers of the angular momentum. The Yi m 's refer to spherical har- 
monics and S R is the MT radius. The numerical orbital (f>^(E,r) is augmented 
inside the sphere by a renormalized spherical Bessel function, 

(2Z + 1)!! 

?*) 

whereas outside the spheres, a renormalized spherical Hankel function is added: 

H kA (r)=i l Y l J-!0^ h l (kr). (3) 

Here, hi = ji — ini is a linear combination of spherical Bessel and Neumann 
functions. Thus, the tail of the MT orbitals at r > Sr is the solution of the 
Hclmholtz equation with zero kinetic energy The potential parameters p A are 
chosen so as to ensure that the wave function be continuous and differcntiable 
at the sphere boundary. 

In this approximation the potential V xc is either von Barth-Hcdin (vBH)[13] 
or Vosko-Ceperly-Aldcr [14] type at the Local Spin Density Approximation 
(LSDA) level, whereas at the Generalized Gradient Spin Approximation (GGS) 
level— which is a functional containing a density gradient correction— the poten- 
tials used the Langreth-Mehl-Hu [15] or Perdew-Wang type [16]. The exchange- 
correlation potential used in the present calculations corresponds to that con- 
taining the von Barth-Hedin parametrization, whose general form is 

v ° c = A(r s )(^Y +B(r s ), (4) 



JkA{r)=i l Y lm \s J" ji(kr), (2) 
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where A(r s ) and B(r s ) are analytical functions. 

2.2 Geometry Optimization 

The ternary system unit cell was optimized by means of ab initio methods, based 
on a DFT treatment within the Local Spin Density Approximation (LSDA) and 
with the Generalized Gradient Spin (GGS) approximation. In the CASTEP 
program [10] the wave function is expanded in planes waves for the valence 
electrons, whereas the core electrons (bound to the nucleus) are taken into 
account by means of their effective interaction on the valence electrons, in the 
form of pseudopotentials which are added to the corresponding Kohn-Sham 
Hamiltonian. 

The pseudopotentials used in this work were generated by Vanderbilt [17] in 
the Klcinman-Bylander [18] representation. The parametrizations of polarized 
spin developed by Perdew-Zunger [19] and Perdew-Burke-Ernzerhof [20] for the 
exchange and correlation energy were used. The conjugated gradient method is 
employed to relax nuclear positions. The sampling of the Brillouin zone was of 
7x7x8 and 9 x 9 x 11 for the LSDA and GGS approximations, respectively, 
using the scheme of Monkhorst-Pack [21]. The cutoff energy for the plane waves 
was of 400 eV approximately. Self-consistency in the calculations was attained 
whenever the total energy changes were < 5 meV, which corresponds to a cri- 
terion of reasonably good convergence. 

To perform geometry optimization we let the lattice constants a, b and c vary as 
free parameters— though they are expected to undergo changes no larger than 
5% relative to the experimental values. The crystal energy is minimized with 
respect to the degrees of freedom by taking into account the calculation of 
Hellmann-Feynman forces in the atoms and the components of the stress tensor 
[22]. Finally, the utilized optimization criteria were 0.00002 eV, 0.0010 A and 
0.050 eV/A for energy change, quadratic mean displacement, and quadratic 
mean force per atom, respectively. 

2.3 Mathematical Structure of the approximation 
LDA+U 57C 

In the LSDA+U [9] method the electrons are separated into two subsystems 
(i and ii). For Ce, (i) delocalized s,p and d electrons, which are described 
by an orbital- independent one-electron potential (the LSDA potential), and (ii) 
localized 4/ electrons, for which we take into account the orbital degeneracy and 
a Coulomb interaction of the form \U^ Jl7 _ tlJ i n a n a i , where n a is the /-orbital 
(or d-orbital) occupancy. For Fe, (i) delocalized s,p electrons, and (ii) localized 
3d electrons. 
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The Hamiltonian for the spin orbitally degenerate systems is of the form 

H = ^ima^im'a + ^ ^ ^ ^ "imj"im'<T 

z, j m.ra' cr i m^m' cr 

+ "2" ^ ' ^ ' "imj^im 1 -cr (5) 

where Cj m 'a is an electron annihilation operator with an orbital index m and 
spin ct(= a,/3). In the lattice site i,tY- m are the hopping integrals and fiim'-a 
is the number operator of the / (or d) electron at site i, orbital m with spin 
a. The first term in Eq. (5) describes the hopping of electrons between lattice 
sites i and j; the interactions between the localized electrons are described by 
the second and third terms, where U and J represent the on-site Coulomb and 
exchange interaction, respectively. 

If we want to correct the LSDA functional for localized electrons we must first 
extract their LSDA treatment to avoid double count of the interaction. The spin 
density functional theory assumes a local exchange-correlation potential which 
is a function of the local charge and spin densities, so, fluctuations around the 
average occupations are neglected. In the mean field approximation (MFA) we 
can write 

where n ma is the mean value of h ma and n a — J2 m nma ■ By introducing this 
approximation in Eq. (5), we obtain the expression for potential energy in the 
mean field approximation, 

E MF = — J ^] y~] n ia (n ia - n vma ) + y n iv n i-cr ( 7 ) 



Solovycv et al. [23] propose to extract an energy function from the total number 
of electrons per spin which would act as the LSDA potential. Such expres- 
sion can be obtained from Eq.(5) in the atomic limit where occupation of the 
individual particle m„ is either or 1: 

id icr 

This energy is now subtracted from E MF to obtain that associated with the 
total energy for localized states: 

i<7 

ima 
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The fraction of the potential acting on the localized orbital (ma) is found by 
differentiating Eq. (9) with respect to the occupation number ni ma : 



dAE 



= AV, 



(U -J){\- »*. 




(10) 



An orbital dependent one-electron potential is thus obtained. 



3 Results and discussion. 



3.1 First Stage: LMTO-ASA calculations 



Results obtained in the optimization of the intcrmctallic compound CeFeGe3 
are shown in Table II. The LSD and GGS approximations yield values below 
the experimental parameters, which may be due to the fact that neither ap- 
proximation removes the error introduced by the self-energy interaction term, 
arising from double count in the Hamiltonian. Furthermore, this material does 
not present a trend similar to that of semiconductors since otherwise, the LSDA 
and GGA results would give values remaining below and above the experimen- 
tal parameters, respectively. A similar situation is found in the equilibrium 
properties of the plutonium (phase <5-Pu)[24], in which, results obtained with 
LDA and GGS give numbers that remain below the corresponding experimental 
values. We work with the CASTEP program at GGS level employing 446 K- 
points in the Brillouin zone. On the other hand, the density of states (total and 
partial) and band structure were obtained via the LMTO-ASA methodology 
within the DFT theory and the LSDA approximation, using the von Barth- 
Hedin parametrization at the optimized unit cell geometry. The corresponding 
results are compiled in Table I. 

The total density of states (DOS) is given in Fig. 1, whereas the partial DOS 
associated with cerium /-states and iron d-states is illustrated in Figs. 2 and 
3, respectively. They all indicate a metallic behavior for the material here 
analyzed, since we have the Fermi level slightly displaced from the central band. 
This is also supported by looking at the plot of band states (BS), as shown in 
Fig. 4, where a dense concentration of bands due to the cerium / states occurs 
around the Fermi level. Furthermore, the fact of having almost horizontal bands 
points to the characteristic behavior of a heavy fermion compound (in our case, 
7 = 150 mJ-molK 2 ) [1,5]. 

In the literature it is reported that the greater magnetic contribution stems 
from cerium (being the magnetic impurity responsible for the Kondo effect); a 
fact that is corroborated in the present theoretical calculations (see Table I). 
Once optimized the geometrical parameters, the material here analyzed turns 
out to be one of those classified as an intermediate valence compound. In our 
calculation, Ce has a magnetic moment of —0.00132 n.p whereas the compound's 
total magnetic moment is —0.0010084 up. Therefore, according to the criterion 
proposed by Vildosola and Llois [25], i.e., by using calculations at the LSDA 
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level and the exchange-correlation potential of Perdew-Wang, they propose that 
the materials can be classified as follows: 



Itinerant 


if 


MCe 


- o, 


Intermediate valence 


if 


MCe 




Magnetic 


if 


MCe 





An intermediate valence system is usually defined in the literature as one having 
a noninteger average number of / electrons. Note that the charge on the Ce atom 
according with the LSDA calculation is of —2.4654 a.u. with a magnetic moment 
of —0.00132. This means that the normal atom configuration of [ATe]4/ 1 5d 1 6s 2 
is being added with almost 3 electrons whose alpha and beta spins pair off one 
another, which results in an overall moment of nearly zero magnitude. 



3.2 Second Stage: LSDA+U calculations 



As can be inferred from the analysis presented in the previous stage, the 
DFT theory— because of its intrinsic formulation— it cannot deal properly 
with strongly correlated systems, therefore, we have chosen to resort to the 
LSDA+U 5/C approximation, developed by Anisimov et al. [9], to calculate the 
electronic structure of the intermetallic compound CeFeGe 3 , as implemented 
in the Stuttgart's TB-LMTO-ASA code version 47. For this calculation, we 
employed the parameters reported in the literature: U = 5.4 eV [11] for cerium 
and U = 2.3 cV (J = 0.9 cV) [12] for iron. 

The resulting density of states in this case is plotted in Fig. 5, whereas the 
partial DOS for cerium 4/-states and iron 3d-states are depicted in Figs. 6 
and 7, respectively. By virtue of the approximation used in this stage, the 
latter include effects arising from strong correlation, not accounted for in those 
obtained by means of a conventional DFT calculation. The corresponding band 
structure is displayed in Fig. 8. 

In the display of total DOS the greatest contribution comes from the cerium 
4/-states where an increased splitting within the energy bands— in addition to a 
relative overall shift around the Fermi level— can be seen in Fig. 5. This feature 
also shows up in the partial density of states, where a large band energy splitting 
arising from the cerium 4/-states can be appreciated in Fig. 6 (being 5 eV wide, 
approximately), as compared with the corresponding energy band display (Fig. 
2) obtained in the previous stage, in which, the presence of a nearly single 
peak dominates, located very close and above the Fermi level. The magnetic 
contributions to this material come partly from both cerium and iron, even prior 
to the introduction of the exchange parameter (previous stage), whereas the 
magnetic moment for germanium is almost null, as seen in Table I. Furthermore, 
the material presents metallic behavior. When carrying out calculations with 
an exchange parameter J equal to 0.90 eV for iron there occurs a reduction 
of the magnetic moment in a 42% rate, which probably follows from a slight 
localization effect of electrons on the iron 3d-shell. On the other hand, the high 
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concentration of energy bands occurring around the Fermi level is consistent 
with a metallic behavior. In fact, a greater density of bands is observed in 
this stage (see Fig. 8) as compared to that observed without introducing the 
Coulomb parameter (see Fig. 4). Such feature— in this stage— points to a typical 
characteristic of heavy fermion materials. 

According to the criterion mentioned in the first stage, proposed by Vildosola 
and Llois [25], and extended to handle the total magnetic moment per unit 
cell, one can notice that the magnetic moments for cerium, — 1.1407/i^j, iron, 
(0.6147/i^), and —0.00055^ f° r the empty spheres, give altogether a total mo- 
ment of — 0.552/1,3. This falls within the classification corresponding to an mag- 
netic material. Also, the charge distribution shown in Table I (empty spheres 
introduced by the ASA condition), points to a possible covalent character among 
the atoms of cerium, iron and germanium, although this manifests iself very 
weakly, despite the fact of accounting for correlation effects in the theoretical 
calculation. Comparing the charge of the Ce atom in the LSDA(vBH) and 
LSDA(vBH)+U cases, it is seen to be similar for both: —2.4654 a.u. and 
—2.3847 a.u., respectively. However, the U interaction gives rise to a complete 
change in the spin behavior: whereas in the absence of U almost 3 electrons 
are added to the neutral atom— although the overall magnetic moment is nearly 
zero— now the salient effect in the presence of the U interaction is to align the 
electrons, yielding an effective magnetic moment whose magnitude is — 1.1407/1,3. 

Having no charge difference in the two cases, where the U interaction acts only 
on the / electrons, the filling of the /-shell proceeds via splitting of the a and (3 
levels by the energy U and the occurrence of electron alignment on the /-shell 
only. The remaining charge distributes over the other shells: s, p and d. 

4 Conclusions 

The calculations performed by means of the LMTO-ASA approximation, within 
the DFT theory, lead to results that are similar to those reported in the liter- 
ature, in particular, the obtained magnetic contribution of the compound here 
analyzed, which practically corresponds to that of a nonmagnetic material. On 
the other hand, the partial DOS, together with the band structure, show a 
metallic behavior for the latter. Furthermore, results obtained with calculations 
carried out where the Coulomb parameters U and J (for iron) arc introduced, 
also favor a metallic behavior and, in addition, a heavy fermion character for 
this material. The two-stage analysis performed in the present study also in- 
dicates a small charge covalent character. The pronounced magnetic moment 
reduction occurring in iron is here ascribed to an electronic cloud localization 
on the 3rf-shell of this atom, which arises as a direct consequence of taking into 
account strong electronic correlation effects. 
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Table and figure captions 



Table I: Parameters utilized in the electronic structure calculations of the 
ternary compound CeFeGe3 and theoretical data obtained by using the TB- 
LMTO-ASA approach. 

Table II: Cell parameters optimized by means of the CASTEP program. 

Fig. 1: Total density of states of CeFeGe3 obtained by the LSDA-vBH approx- 
imation. 

Fig. 2: Partial density of states (Ce 4/-states) obtained by the LSDA-vBH 
approximation. 

Fig. 3: Partial density of states (Fe 3d-states) obtained by the LSDA-vBH 
approximation. 

Fig. 4: Band structure obtained by the LSDA-vBH approximation. 

Fig. 5: Total density of states obtained by the LSDA(vBH)+U approximation 
using the parameter U for cerium and parameters U and J (exchange) for iron. 

Fig. 6: Partial density of states (Ce 4/-states) obtained by the LSDA(vBH)+U 
approximation. 

Fig. 7: Partial density of states (Fe 3d-states) obtained by the LSDA(vBH)+U 
approximation. 

Fig. 8: Band structure obtained by the LSDA(vBH)+U approximation using 
the parameter U for cerium and parameters U and J (exchange) for iron. 
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Table I 



Atoms in the 


Crystallographic 


MT sphere 


Charge" 


Magnetic" 


Charge 6 


Magnetic 6 


CeFeGe 3 


Positions 


radii 


(a.u.) 


Moment 


(a.u.) 


Moment 






(a.u.) 




(Ms) 




(Ms) 


Ce 


(0.0, 0.0, 0.0) 


4.1077 


-2.4654 


-0.00132 


-2.3847 


-1.1407 


Fe 


(1.0,1.0,0.66) 


2.4566 


-0.3027 


0.00033 


0.3339 


0.6147 


Gei 


(0.5,0.0,0.25) 


2.5320 


1.1212 


-0.000005 


1.1208 


-0.01274 


Ge 2 


(1.0,1.0,0.42) 


2.6284 


1.064 


-0.000013 


1.0499 


-0.02511 


Ei* 


(0.1224,0.1224,0.5684) 


1.1281 


-0.1396 


-0.000001 


-0.1432 


-0.00055 


E 2 * 


(0.1224,-0.1224,0.5684) 


1.1281 


-0.1396 


-0.000001 


-0.1432 


-0.00055 


E 3 * 


(-0.1224,0.1224,0.5684) 


1.1281 


-0.1396 


-0.000001 


-0.1432 


-0.00055 


E 4 * 


(-0.1224,-0.1224,0.5684) 


1.1281 


-0.1396 


-0.000001 


-0.1432 


-0.00055 



* Positions of the empty spheres that fulfill the ASA condition within the LMTO-ASA 
approximation. 

"As obtained for the parameters optimized by means of the LSDA-vBH approximation. 
6 As obtanied via calculations performed with the LSDA(vBH)+U aproximation. 



Table II 



Experimental cell 


Optimized cell 


% of error 


Optimized cell 


% of error 


parameters 


parameters 


Exp. vs LSDA 


parameters 


Exp. vs GGS 


(A) 


(A) LSDA 




(A) GGS 




a=b=4.332 


a=b=4.1767 


3.71 


a=b=4.234 


2.314 


c=9.955 


c=9.5981 


3.72 


c=9.73 


2.312 
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